Mon. Not. R. Astron. Soc. 000, 000-000 (2008) Printed 15 April 2008 (MN MfcjX style file v2.2) 



Oscillating pulsar polar gaps 



oo 
O 

o 

S-l 

Oh 
< 



Oh 
6 



> 

o 
o 

(N 

o 

oo 
O 



X 



Qinghuan Luo and Don Melrose 

School of Physics, The University of Sydney, NSW 2006, Australia 



Received in original form September, 2007 



ABSTRACT 

An analytical model for oscillating pair creation above the pulsar polar cap is pre- 
sented in which the parallel electric field is treated as a large amplitude, supcrluminal, 
electrostatic wave. An exact formalism for such wave is derived in one-dimension and 
applied to both the low-density regime in which the pair plasma density is much lower 
than the corotating charge density and the high-density regime in which the pair 
plasma density is much higher than the corotating charge density. In the low-density 
regime, which is relevant during the phase leading to a pair cascade, a parallel electric 
field develops resulting in rapid acceleration of particles. The rapid acceleration leads 
to bursts of pair production and the system switches to the oscillatory phase, corre- 
sponding to the high density regime, in which pairs oscillate with net drift motion in 
the direction of wave propagation. Oscillating pairs lead to a current that oscillates 
with large amplitude about the Goldreich- Julian current. The drift motion can be 
highly relativistic if the phase speed of large amplitude waves is moderately higher 
than the speed of light. Thus, the model predicts a relativistic outflow of pairs, a 
feature that is required for avoiding overheating of the pulsar polar cap and is also 
needed for the pulsar wind. 
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1 INTRODUCTION 

One of the central problems in pulsar electrodynamics is the 
production of the relativistic electron-positron pair plasma 
in whi ch coherent rad io emission is thought to be pro- 
duced |Sturrockl Il97ll ). It is widely believed that parti- 
cles are accelerated to ultra high energy along open field 
lines, leading to a cascade producing the pair plasma. The 
pair cascades should produce detectable X-rays and gamma- 
rays in the case of fast rotating, young pulsars and mil- 
lisecond pulsars (|Thompsonl l200ll h Different acceleration 
locations in the pulsar magnetosphere have been postu- 
lated, with a common feature that the acceleration re- 
sults from a large scale electric field on open field lines 
that extend beyond the light cylinder. Acceleration regions 
in a pulsar magnetosphere are referred to as 'gaps'. A 
widely-discussed acceleration region is near the polar cap, 
for which there is a class of acce l eration models called 
polar gap models dSturrockl Il97l1 ; lArons fc Scharlemannl 
1 19791 ; lHarding fc Muslimov 199a ) ■ Particle acceleration near 
the polar cap is of particular relevance for pulsar ra- 
dio emission since the radio observations suggest that for 
many pulsars the emission originates deep inside the pulsar 
magnetospheres dBlaskiewicz, Cordes fc Wassermanl Il99ll ; 
lEverett fc Weisberd [200^ 7 Other acceleration sites include 
regions in the outer magnetosphere near the null surface, 



referred to a s outer gaps (Cheng. Ho fe Rudermanl Il986l ; 
lRomani|[l996l : iHirotanil 12006) and a variant, referred to as 
slot gaps, charact erized by a l o ng, thin region along the las t 
open field lines (|Aronsl 1 19831 ; lHarding fc Muslimovl 120051 ). 
Here we concentrate on the polar cap region. 

The conventional polar gap models were developed on 
the basis that the system can settle into a steady state, such 
that all physical quantities in the pu lsar's corotating frame 
can be regarded as time- independe nt l|Arons fc Scharlemannl 
1 19791 ; lHarding fc Muslimovl Il998ft . This time-independent 
assumption is not realistic in practice because it ignores in- 
ductive electric fields. There is a strong argument that a 
time- dependent inductive field plays a central role in pulsar 
electrodynamics. A global curren t must be present in the 
system to form a current closure fS charlem ann fc Wagoner! 
1 19731 ; iMichelll 19751 ; ICheng fc Rudermanlll976l : IShibatalll99lD 
and the steady-state assumption requires that the current 
density, J, balance V x B/no- Moreover, the global electro- 
dynamics requires that the current be determined globally, 
rather than by local p rocesses near the polar cap l|Shibatal 
Il99ll ; iTimokhinl l200d l. The global and local requirements 
on the current density are generally incompatible without 
an additional source of charge and current, and any result- 
ing mismatch between J and V x B/flo implies a time- 
dependent electric field. Such a mismatch and associated 
inductive electric field seem unavoidable. 
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Although the idea that a cascade above the polar cap 
is intrinsica l ly tim e-dependent w as suggested much earlier 
bvlSturrockl (|l97ll ) and adopted in lRuderman fc Sutherland! 

spark model, there are no quan titative models 
that t ake into account the time dependence. iLevinson et al.l 
(2005) recently discussed an oscillatory gap model to illus- 
trate the time dependent nature of the polar gap in the one- 
dimensional approximation. In this model, induction cur- 
rents due to temporal changes in the system are included 
and particle acceleration is shown to settle into an oscilla- 
tory state, similar to a large amplitude wave. The model 
has several limitations. First, escape of particles is included 
only implicitly, and escape needs to be made explicit to be 
consistent with observations of pulsar winds. Second, the os- 
cillatory model is based on numerical integration of the rel- 
evant fluid equations together with Maxwell's equations in 
one dimension and due to the limit of numerical calculation 
the oscillations could be followed only for a limited number 
of periods. Third, the assumption that the oscillations are 
purely temporal is unrealistic, and needs to be generalized 
to allow outward propagating waves. 

In this paper we adopt a different approach, treating the 
oscillations as a large-amplitude electrostatic wave (LAEW) 
in a cold pair plasm a. Although LAEW in a cold electron 
gas was discussed in lAkhiezer et all l|l975h . to our knowl- 
edge there has been no discussion of the case of a strongly 
magnetized, electron-positron pair plasma. Relativistic mo- 
tion of a single particle in a LAEW was d iscussed and ap- 
plied to pulsar emission by iRowd (Il992al lbh. However, in 
the single-particle treatment, the wave was assumed to pre- 
exist and feedback of the particle's motion on the wave field 
was not taken into account. Here we treat the electrons 
and positrons as cold fluids and determine particle accel- 
eration by solving simultaneously both fluid equations and 
Maxwell's equations without making an a priori assumption 
of time-independence. We include pair creation in our equa- 
tions, but neglect it in deriving analytic solutions. 

In Sec [2] we outline the fluid formalism for time- 
dependent electrodynamics that includes inductive electric 
fields. Analytical solutions for LAEW are described in Sec[3] 
Inclusion of pair production in LAEW is discussed SecUand 
the low-density limit is discussed in Sec [5] 



2 TIME-DEPENDENT FORMALISM 
2.1 Fluid equations 

We consider a two-component cold fluid consisting of elec- 
trons and positrons; the fluid number density and velocity 
are denoted by N± and v±, where the subscripts ± corre- 
spond respectively to the positron and electron components. 
In the observer's inertial frame, a pulsar rotates with an an- 
gular velocity f2 = 2tt/P, where P is the pulsar period. Pro- 
vided that sufficient charge density is available, the corotat- 
ing electric field is set up so that charged particles corotate 
with the star. This electric field can be eliminated by choos- 
ing a frame corotating with the star. Well inside the light 
cylinder, the corotating frame and observer's inertial frame 
are connected by a local Galilean transformation with veloc- 
ity wj! = S~l x r , where r is the radial vector directed from the 
star's center to a point of interest. One has |«h/c| ~ r /Rlc, 



with Rlc — c/fi the light-cylinder radius. Thus, the effect 
on the rotation of the fluids can be neglected if the region 
concerned is close to the polar cap where r / Rlc <C 1. 

The relevant fluid equations, the continuity equation 
and equation of motion, can be written down as 

dN ± 



dt 
d_ 
dt 



+ V • (N±v±) 



v± ■ V I u 



Q 

2N± 



u±, 



(1) 



(2) 



where q ± is radiation drag, Q is a source function due to 
pair production and u± = j±v±/c is the particle's dimen- 
sionless momentum. The Lorentz force is absent because all 
the particles are assumed to be in the ground Landau state. 
The current and charge densities are 



J = e >J sv a N s , p — e >J sN s 



(3) 



with s = ±. The fluid equations are supplemented by 
Maxwell's equations, written in the corotating frame. The 
two relevant Maxwe ll equations involve the charge dens ity 
and current density (jFawlev. Arons fc Scharlemannlll977l ) 



V ■ E = 



£"() 



(P - 



V x B = p, (J - J R ) + — — , 

" 2 at 



(4) 
(5) 



where p GJ = e [-2il- B + x r)- (V x B)} is the Goldreich- 
Julian (GJ) density. Here we ignore the general relativistic 
effects such as frame dragging (cf. Sec 5.2). In ([5j) Jr is 
a vectorial sum of all the remaining terms that are small 
fo r t/Rlc <S 1; the full expression for J r is given by (A5) 
in lFawlev. Arons fc Scharlcm anrJ l|l977h and it is neglected 
here. Equation @ describes a noncorotating electric field 
arising from deviation of charge density from the GJ den- 
sity. In steady-state models, equation (JSJ) is (implicitly) as- 
sumed to be satisfied trivially. However, it plays a central 
role here in determining the inductive field arising from a 
current mismatch. 

The inclusion of the inductive field distinguishes the 
model considered here from steady-state polar-gap mod- 
els in which only Poisson's equation (i.e., Eq fij) is rel- 
evant and the paralle l elect r ic field is treated as stati c 
l|Arons fc Scharlemannl 1 19791 ; lHarding fc Muslimovl 1 19981 ). 
The static assumption is incompatible, in general, with 
the constraint imposed by a global current. It has long 
been recognized that circulation of a global current 
plays a crit ical role in dissipation of rot a tional e nergy 
of pulsars (IScharlemann fc Wagoner] Il973l ; iMichell Il975l ; 
ICheng fc Rudermanlll976l ). Such current, denoted by Jo, 
should be determined globally. The simplest case is where 
the global current is assumed to be a constant. Assuming 
Jo|| = (V x -B)||, where || represents projection along the 
magnetic field, the steady-state assumption implies that the 
local current exactly matches the global current, Jn = J || 
(when Jr is ignored), everywhere. This assumption is not 
realistic, for exampl e, due to pair creation changing the cur- 
rent density locally l|Levinson et al.|[2005l ) . In our oscillatory 
model, the parallel electric field is treated predominantly as 
an inductive field due to oscillation of the current about the 
global constant direct current. Relevant solutions are oscil- 
latory and should be treated as large amplitude waves. 
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2.2 Large amplitude waves 

We outline a general approach for deriving a time-dependent 
solution in which oscillations are considered as a large am- 
plitude wave propagating in the direction k at a constant 
phase speed /3v (in units of c). We assume that oscillatory 
quantities are functions of 



X = w, 



w (ft 



(6) 



where w GJ = (e 2 NGj/eoTrie.) 1 is the plasma frequency 
at the GJ number density, Nqj = 2eo^-Bo/e, and Bo is 
the surface magnetic field. We normalize the electric field 
as E\\ = eE\\/m e cu) G j , density as N± = N±/Ngj, the 
charge density as r\ = p/ecNcj and the current density as 
j = J/ecNGj- Using d/dt = UGjPvd/dx and V = d/dr = 
~(ujGjn/c)d/dx, one may write fTJ, © and © in the di- 
mensionless forms 



_d_ 

dx 



Q3v-k-P ± )N± 



du± 



±E - 



1± 



2N± 



u±, 



(1 



■ft) 



dE 



(7) 



(8) 



(9) 



where /3 ± = 



9± = q±/(m e c u GJ ), Q = 
y and j = V x B /(cNgjc) which is 



»±/c 

Q/{N G ju) gj ), [Tjij = <5 ; 
assumed to be a constant vector. Poisson's equation (4) and 
the induction equation (5) may be replaced by the equation 
of charge continuity and the induction equation. With the 
charge and current densities functions only of x, the equa- 
tion of charge continuity gives 



0. 



(10) 



Note that in our model both Poisson's equation and the 
induction equation contribute. The steady state models cor- 
respond to the limit j3v —> 0, and hence K • j = const., 
when the induction equation does not contribute. It is only 
in the opposite limit, /3y — * oo, where the oscillations are 
purely temporal, that Poisson's equation makes no contri- 
bution (cf. Sec 2.3). As we are primarily interested in os- 
cillatiory solutions, in deriving these equations we ignore 
inhomogeneity, notably in the magnetic field. This neglect 
is justified provided that the spatial scale for the inhomo- 
geneity is much larger than the oscillation length. The scale 
of the inhomogeneity is of order the radius of curvature, 
which at a height r — R with R — 10 4 m the star's radius, is 
larger than but of order (tRlc) 1 ^ 2 - One has ~ (tRlc) 1 ^ 2 ~ 
2.2 x 10 s {r/R) 1/2 P Q 1/2 m, where F .i = P/0.1 s. For plausible 
parameters, this is much larger than the oscillation length, 
which is of order 



A = 



UGJ \N ) 



0.9f3 v Bl /2 P 



1/2 



f jlV /2 (JL 

UO 6 / \ 10 2 



-1/2 



(11) 



with B s = B /10 s T, N = max{iV + , iV_ } and where we ass- 
sume an oscillation frequency ~ tu v /'y 1 ^ 2 (see Ea l33p . For the 
numerical example in , the condition is satisfied provided 
that (3y < 2.4 x 10 . We are only concerned with the open 



field line region which can be regarded as a flux tube with a 
conducting surface, defined by the last closed field lines, that 
separates the region from the closed field line region. Eq 
((9| are then valid only when the wavelength is much shorter 
than the transverse size ~ (r / 'Rlc) 1 ^ 2 t ^ 258P ~ 1 1//2 m. 
Such short wavelength approximation implies that the limit 
0v — > oo is not applicable for pulsars. Nonetheless, such 
limit is also discussed here as it simplies the formalism from 
which the basic properties of LAEWs can be derived and 
compared to a more general case where f3v is finite. 

Integration of Eq (Q) yields an exact form for the density 

(0v ~ k ■ f3 ±0 )n ± +F Q 1 fx * , , 



N ± = 



Pv - K- (3, 



where n± — N±(0) is the initial density at x = 0, /3 ±0 are 
the initial velocities (in units of c), and Fq is a cumula- 
tive flux arising from pair creation. A wave is classified as 
superluminal if /3v > 1, luminal if f3v = 1 and subluninal 
if Pv < 1. We do not consider the subluminal case here. 
For superluminal and luminal waves, N± is always positive. 
The number density remains approximately constant in a 
superluminal wave in the limit /3y — > oo. 

2.3 Current-charge invariant 

The charge continuity equation (|10p implies an invariant k ■ 
j — /3yr) — const. Denoting the dimensionless current density 
by i|| = K ' ii this conservation law implies 



h(x) - Pwix) = j\\(Q) - Pvv(o), 



(13) 



where \ = corresponds to the initial conditions, with 
J||(0) = /3+on+ — /3_on_, rj(0) = n + — n_ . The steady state 
condition corresponds to the special limit j3y — > 0, and in 
this limit the induction equation is satisfied trivially with 
ill = Joil • The constant current j is interpreted as the global 
direct current, assumed to be determined by global condi- 
tions and to be a free parameter in the theory. In steady 
state models, only Poisson's equation is relevant and the 
constant current j u does not appear explicitly. For j3y ^ 0, 
it is convenient to write rj GJ — p CJ /eNGj, for the sign of 
the GJ charge density, and to write the right hand side of 
Eq (HJ) in the form 



i||(0)-/3vr?(0)=j || -PvVc 



(14) 



The opposite limit /3y — > oo corresponds to purely temporal 
oscillations. In this case (|14|l requires that the charge density 
equal the GJ density, so that Poisson's equation is satisfied 
trivially. We are interested in the general case < f3v < oo. 

Assuming a strong magnetic field approximation, a 
plausible simplifying assumption is that n is directed along 
the magnetic field. The problem then becomes one dimen- 
sional, involving only projections of the relevant equations 
along the magnetic field. Of particular significance is j j| > 
which would be identically zero if the magnetic field were 
dipolar. The global requirement for j n 7^ implies a 
nonzero azimuthal magnetic field B<$, ~ (r / Rlc)jo\\B <g; B 
(see further discussion in Sec 3.4). Although j || cannot be 
determined locally, it is plausible to assume that it has the 
same sign as tjgj, i.e., j ii < for CI ■ B > and j || > for 

n ■ b < 0. 

One may use (|13[) and (|14[) to express Pv in terms of 
the initial density and velocity: 
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Pv 



Jo II 



/3+on + + P-on- 



— 71+ + 71- 



(15) 



provided that rj(0) — n+ — n_ 7^ i/gj- The following three 
cases are of interest: the initial charge density matches the 
GJ density 8-q = t)gj — T](Q) = 0, the initial charge density 
has a small deviation from the GJ density \Srj\ <C 1, and 
pair density is much lower than the GJ density n± <C 1. In 
the first case, one must have j\\ (0) = jo\\ and thus, the phase 
speed is not constrained by flD - Generally, the second ap- 
plies to polar cap regions where both frame-dragging and 
field line curvature may cause a smal l deviation of a charge 
density from the corotation density IIArons fc Scharlemannl 
Il979l : IShibatalll997l : lHarding fc Muslimovll 19981 ). Eg fify inT 
plies \f} v \ > 1 for \6r)\ < |1— A>|, jo\\ = r\ GJ and P = /3 +0 w 
P-o- Therefore, the LAEW considered here should be super- 
luminal. In the third case, one has 0v ~ jow/WGJ - , the wave 
is superluminal for j « > t]gj and subluminal for j u < -qcj- 



3 LARGE AMPLITUDE, ELECTROSTATIC 
WAVES 

We consider large amplitude, electrostatic waves in a high 
density regime where the density of pairs is much higher 
than the GJ density, n± S> 1, and pair creation is absent. 
This regime is especially applicable when the system under- 
goes a brief burst of pair production leading to a pair plasma 
with n± ^> 1 and sets up oscillations. 

To concentrate on the basic physics of such large am- 
plitude waves we ignore the radiation drag, q± = 0. Eqs ((9} 
and © reduce to the following simple forms: 



(3 s o)n s 



d X PV Pv ^ PV~ P: 



d 



PvPi 



= TE|| 



d X {i -Pi) 1 ' 2 

Integration of ()17[) leads to the following invariant: 
1-PvPs 1-PvPso 



Sri 



(l-A 2 ) 1 / 2 



Sri 



(i - P 2 ) 1/2 



6 



(16) 



(17) 



(18) 



One needs to consider only one component of the fluid, say 
the electron (— ) and the solution for the other component 
can be derived using (|18|) . Plots of P+ as a function of P- 
are shown in figure [TJ for luminal and superluminal waves. 
Using the notation 

1 - PvP- 



£ 



(l-/3i)l/2' 

the velocity can be expressed in terms of £: 

Pv-i(p 2 v+e-i) 112 



P-(0 = 



Pi+e 



(19) 

(20) 
(21) 



P+(0 = M&-0- 

Notice that the two velocities (|20p and ()2ip are related by a 
transform £ — > £1 — £. Exact solutions to (|16p and (|17p are 



/ 



$l/2(£) 



= X, 



1/2 



(22) 



(23) 



with 

*(£) = ^o + #- 

Pv 



(£ - £o)jo - (Pv - P+oMOn+ 



-08v-/3_o)(7-(O-7-o)n- 



(24) 



9 = 



/?+(£') 



, o Pv-P + {t') 
1 ' 



PI -I 



£0 - £ - /3v 



(/# + (6 - £) 2 



1/2 



•(/3v + (a-eo) 2 



1/2' 



(25) 



where £0 = £(0) and 7-0 = 7-(£o) is the initial Lorentz fac- 
tor of electrons. The initial electric field is Eo = ±$ 1//2 (£i) = 
d£,/dx at x = 0, where the sign is determined by the sign of 
d£,/dx at x = 0. Although we are interested in superluminal 
waves, the calculation up to this stage applies to any j3y in- 
cluding the special case j3v = 1. For f3v ~ 1, (|25p simplifies 
to 9 = [(£1 - £) _1 ~ (6 ~ £o) _1 + £0 - £]/2. The function 
<?(£) is shown in figure [3] for superluminal waves f3v > 1. 

3.1 Condtions for an oscillatory solution 

The condition for existence of an oscillatory solution for (|16p 
and Jl7|) is $(£) > 0. For n± > 1, (JUJ simplifies to 



$ « £ 2 - (2/(3v)(Pv - P-o)(g + 1-- 7o)n- > 0, 



(26) 



where (|15p is used to eliminate n+. An appropriate approx- 
imation to g(£) can be derived as follows. In each oscilla- 
tion a particle's velocity is highly relativistic except for a 
short phase when it is briefly nonrelativistic and changes 
sign. In this case, (|25p can be approximated by two straight 
lines that cross at £ w £1. For /3±o — (solid and dashed) 
the minimum is located at £ = £1 = 2. Specifically, for 
(£1 - £o) 2 » Pv ~ 1 and (£1 - £f » - 1, °ne has the 
following approximation 

£0 - £1 



+ 



£-£1 



£>£i, 



£ < £1, 



1 + M£i - £o)/V Pv + 1 

with h(x) = 1 for 1 ^ and h(x) = —1 for x < 0. In 
the large /3y limit g(£) is symmetric about the vertical axis 
at £ = £1. As the phase velocity approaches the luminal 
limit Pv — * 1 the right-hand side of the curve steepens and 
approached the vertical axis. In figure [3] (|27[1 is shown as 
dotted lines, which gives a quite good approximation to the 
exact numerical result (solid lines). 

The condition (|26p leads to upper and lower limits to 
U-, given by 

"max W ^-[(/3v+l)F-£l], 

Ipv 



with 



7-o + 



-2^ 



PvEq 



1)F + £i]. 



£o-£i 



2n-{Pv-P-a) l + h(£i-£ )Pv 



(28) 



(29) 
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P- 




-1 -0.5 0.5 1 



P- 

Figure 1. Plot of /3+ as a function of j3—. Upper: j3v = 1- The 
curves from left to right correspond to fi-j-g = 0, 0.5 and 0.9, 
respectively. In each cases, the physical range corresponds to a 
range from where a maxima occurs to the rightmost. Lower: /3y = 
1.5 and /3±o = (dashed), j3y = 1.5 and /3±o = 0.5 (dash-dotted), 
/3y = 100 and f3 ±0 = (solid), /3 V = 100 and /3 ±0 = 0.5 (dotted). 




-10 -5 5 10 



U_ 

Figure 2. Plot of u+ as a function of u_ for (3q± = (solid), 
(3q± = 0.9 (dashed) and /3q± = —0.9 (dotted). Because of £ ^ 
0, motions of electrons and positrons are not strictly symmetric 
about u+ — » — U-. 




-40 -20 20 40 



I 

Figure 3. Plot of g(£). The three solid V-shaped curves from left 
to right correspond respectively to j3y = 5 with f3±o = 0.9, 0y = 
5 with j3±Q = 0, and py = 1.1 with j3±o = 0. The minima where 
the two lines meet are located at £i ss — 16 for the leftmost plot 
and § = £i = 2 for the other two. The dotted lines are obtained 
using the approximation J27j which consists of two straight lines 
with gradients l/(/3y — 1) and —l/(f3y + 1), respectively. 



Using (|18p the positron momentum u + can be expressed in 
terms of u_ , which leads to the same upper and lower limits 
for positrons as (|28|l. For luminal waves with (3y = 1, one 
has Umax ~ (l + -Eo/ n -)/2 and tt m in ~ —1- In this case oscil- 
lations skew strongly in the direction of wave propagation. 
For fiy — > oo, ([28} reduces to it max = (r + u + o + u-q)/2 
and M m in = — (r — u+o — u_o)/2. When u±o ~ 0, oscilla- 
tions become symmetric with u max = — u m i n = T/2. Since 
itmax 7^ u m i n , oscillating electrons and positrons have a net 
drift velocity ~ (it max + it m in)/2 ~ (r — £i)/f3y. An accurate 
evaluation of the drift velocity is given in Sec. 3.3. 

Numerical solutions to E« = (i.e., <3?(^) = 0) are 
shown as contours in figure [4] We express n+ in terms of 
Pv, Jo || : and r\c.j using (|13p and (| 14[) . In the subfigure on 
the left one assumes luminal waves with (3y = 1. Each pair 
of lines defines an upper limit, Umax, and a lower limit, u m i n , 
to a particle's momentum such that one has $ > for 
Umin < u_ < u m ax- The similar upper and lower limits can 
be obtained for positrons. Particles oscillate with their mo- 
menta confined between these two limits. As u ma x 2> |u m in|, 
the oscillations skew strongly in the wave propagation di- 
rection. The dotted lines correspond to a superluminal wave 
Pv — 2. As Pv increases the system evolves toward sym- 
metric oscillations as shown in the subfigure on the right. A 
nonzero initial velocity (/?o± > 0) also shifts the oscillations 
forward (dash-dotted lines). It should be emphasized that we 
intentionally choose a moderate n_ here to illustrate that 
an oscillatory solution can exist even for a moderate pair 
density. 



3.2 Analytical formalism 

To seek an oscillatory solution one may define the oscillation 
periodicity (in units of 1/uj gj ) 
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Figure 4. Conditions for Eu = 0. (a) A luminal wave /3y = 1, with n_ = 5, Jqm = t\qj = —1. The solid and dashed lines correspond 
to (3±o = and 0.9, respectively. Each pair of lines defines an upper and lower limits to u, between which <I> > 0. For comparison a 
superluminal wave with /3y = 2 and /3q± = is shown as a dot line, (b) Superluminal waves with different phase velocities f3y = 2 
(dotted), 10 3 (solid), 10 (dashed). We assume [3±o = 0. The dash-dotted lines represent f3y = 10 with /3q± = 0.9. 



where the subscripts a, b label respectively the phases at 
which (3- = ,0 m i n and /3_ = /3 ma x, respectively. The os- 
cillation frequency is then given by uj = 2-kuj gj /T, where 
uj G] ~ (1.5 x 10 11 s _1 )P ~ 1 1 _B8. Using the approximations 
(|27p an analytical solution for a LAEW can be derived. For 
Pv ^> 1, using 7-(£) + <?(£) ~ 2|£|//V, one obtains the 
following approximation to (|26l) : 



<!>(^2„ ( | - £j£ ] ( I- 



(31) 



where the upper and lower signs correspond to £ £i and 
£ < £i, respectively. For /3v 3> 1 one may expand (|31[) in 
l//3v and substitute it for (|30[) . obtaining 



T : 



2n-(j3v - P- ) 



1/2 



2go 



(32) 



A second approximation applies for Eg/2n- 3> 7o±, in 
which case one estimates the frequency as 



V2-n 



An 1/2 



where oj d 



(2n. 



(33) 



;j is the plasma frequency of the pair 



plasma. As an example, one has lu w 10 s for n_ 



1CT 



decreases as w ma x increases, which can be understood as an 
increase in the effective mass of electrons or positrons. 

The condition for <3> > is j£| < & ~ — £a ~ PvEo/4n-; 



this gives u n 



E$/4n. 



Note that such symme- 



try in oscillation is a direct consequence of our assumption of 
large Pv and P±o = 0. The solution for < \ < \ T = PvT 
is found to be 



0<X< jXr, 



2n_ V 

k. + ^ : (* ( 1 / 2 -* 1/2 ).!x^x<x I 



(34) 



where <E> is given by pip . $0 
electric field is obtained as 



$(1), $ a , 6 = The 




0<X<X T /2, 
x)> Xt/ 2 ^X<Xt- 



(35) 



The electric field displays a sawtooth wave form, which can 
be understood qualitatively in terms of the extreme rela- 
tivistic limit. In this limit positrons and electrons are ac- 
celerated in opposite directions, giving rise to a current 
~ 2|/3_|n_ 3> I j || I ~ 1. Thus, the electric field is 
E\\ oc ±|j|||//3v ~ ±2n-x/Pv with ~ 1, which repro- 
duces the sawtooth wave form given by (|35|l . An example of a 
numerical integration of (|22[) and (|23[) is shown in figure[5]for 
pairs with an initial, forward velocity. Figure[6]shows oscilla- 
tions for particles with an initial velocity toward the star. Al- 
though our analytical solution is obtained for n_ 2> 1, here 
we again choose a moderate n_ in the numerical calculation 
to show that our oscillatory solution is also valid for n± ~ 1. 
The wave f orm is similar to that predicted from the numer- 
ical model l)Levinson et al.|[2005h . The characteristics of the 
oscillations, including the periodicity and amplitude, are in- 
dependent of the sign of Eq and are not sensitive to the ini- 
tial conditions Po± provided that |lto±| >C Umax- As no radia- 
tive loss, pair production nor wave damping is included, the 
wave amplitude remains constant. In the figure we assume an 
initial electric field much lower than the typical vacuum field 
£max ~ 3 x 10 6 (B/B C ) 1/2 P " 1 1/2 , where B c w 4.4 x TO 9 T. 
In practice, the initial field should be near the pair creation 
threshold. Assuming the pair production threshold to be -yth, 
one has E « 3 x 10 4 (n_/5) 1/2 (7th/10 B ) 1/2 (see discusssion 
in Sec. 5). Since |lt m in| < |M max j, the oscillating particles 
have a net forward flow velocity. 



3.3 Drift motion 

The drift momentum can be obtained by averaging u_ over 
one period (T): 



PvT 



(36) 
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Figure 5. Dimensionless momentum u (upper) and electric field 
-By (lower) as functions of phase \. The dashed line corresponds 
to oscillations of positrons. As there is no dissipation included, the 
amplitudes of the oscillations are determined by \Eq\. Electrons 



oscillate between u„ 



-5.0 and u n 



13.4, with a net drift 



velocity u ps 3.6. Oscillations of positrons (dashed) skew less in 
the positive direction than electrons. We assume Eq = —10, j3y = 
5, f3±o = 0.9, ioy = TjQj = -1, n- = 5. 




25 50 75 100 125 150 175 20 



X 

Figure 6. As in figure \E\but with /3± = -0.99, E = -50. The 
oscillation pattern is similar to figure[5]but with a small backward 
shift. 



velocity f3 D = u/(l + u 2 ) 1 ' 2 . For E$/2 



El 



4u n 



where we expand $ on 1/ f3v <C 1 and change the integration 
variable to d\ = — d£/ ( I >1/ ' 2 (£). Terms of order \j \\ j/n± <g 1 
and \riGj\/ n± ^ 1 are ignored. Then both electrons and 
positrons are dragged along in the wave at the same drift 

3> 7±o one has 

" " 3n-f3 v " Wv ' ^ 
The drift velocity decreases as the wave phase speed in- 
creases. In the limit /3V — ► oo oscillations are purely tempo- 
ral. The upper and lower limits to the particle's momentum 
are u max = (r+u+o+M-o)/2 and M m i n = — (T— u+o— u_o)/2. 
The drift velocity (|36f) reduces to u = (u+o + «-o)/2. 
For u±o — 0, particles oscillate symmetrically between 
^min ~ — iimai ~ E /4n_ and ttmax- It is interesting to 
note that the proportionality l//3y in (137|l is similar to tha t 
predicted from the single-particle treatment (|Rowej[l992bh . 
However, since the single-particle formalism does not in- 
clude the feedback effect of particles on the wave, it pre- 
dicts a low drift velocity /3n ~ l//3v- Our extact treament 
shows that the drift motion can be highly relativistic with 
7i3 = 1/(1 - /3|,) 1/2 ~ \u\ > 1. As a result, the LAEW can 
drive a relativistic outflow of particles even when the parti- 
cles are initially at rest. Thus, the oscillating gap can supply 
relativistic pairs to the pulsar wind. 

An implication of such drift motion is that polar cap 
overheating can be avoided and so, the model satisfies the 
observational constraint on thermal X-rays from the polar 
caps. The observed relatively low fluxes of thermal X-rays 
imply either that the flux of particles that impact on the 
polar cap is much lower than the G,I flux ~ \pgj\c or that 
acceleration of returning particles is insignificant. The latter 
can be ruled out as the returning particles must be subject 
to the same strong accelerating electric field that accelerates 
forward moving particles. In the oscillatory model, since par- 
ticles are dragged forward by the LAEW and escape to in- 
finity, few particles are reflected back to the star; oscillating 
particles do not impact on the polar cap if the oscillating 
region is located a distance > cT/logj ~ A//?y ~ 0.9 m 
for B = 10 8 T and P = 0.1s. So, the model can satisfy 
the observational limit to the thermal X-ray observations 
from t he polar cap. It is worth comm enting that by con- 
strast, iRuderman fc Sutherland! (Il975l )'s 'vacuum sparks' 
model, which is instrinsically time dependent, predicts a 
much larger thermal X-ray flux than the observational limit. 

3.4 Currents 

The current j(x) can be derived using (|34|l . Assuming /3y 3> 
1, we have 



0v 



P-n-) 



(38) 



For half the phase of an oscillation, the electrons and 
positrons are accelerated in opposite directions, and then 
these directions reverse, with only a brief phase in which 
the motion is nonrelativistic. Hence, one has /3+ ~ — /3_ ~ 1 
for nearly all phases. The final term in Eq (|38p is generally 
small, and on neglecting it, the current is j\\ w ±2n_ ex- 
cept for the short phase where it switches sign. Thus, the 
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Figure 7. Current j\\ vs \- The parameters are as in figure[5] The 
square wave form can be understood as that in each oscillation 
particles stay most time in the relativistic regime in which j3± ~ 1. 



oscillating current has a square wave form with an ampli- 
tude Ijmaxl ~ 2n_ 3> |jo|||- A numerical calculation of j\\{x) 
is shown in figure [7] Note that as we consider only electro- 
static waves, so that the oscillating current induces electric 
fields only; there is no oscillating magnetic field. 
The average current is given by 



ill ~ P+n+ - f3-n- w (3(n+ 



(39) 



where j3± is the mean velocity, i.e., drift velocity. The second 
approximation is derived for n± 3> 1 and hence /3+ « /3_. 
The mean charge density is given by fj = rjcj + (j\\ —jo\\)/Pv- 
The system tends to settle into a state where fj — t)gj and 

in = ion • 

3.5 Wave dissipation 

So far we neglect dissipation in obtaining our solution for 
LAEW. This is justified if damping (or growth) is weak in 
the sense that the change in the LAEW in an oscillation 
period can be treated as a perturbation. Damping occurs 
through radiative losses, which include curvature radiation, 
RICS, two-s tream instabili ty, and linear acceleration emis- 
sion (LAE) ljMelroseHl97i l. In the first and second mech- 
anisms, energy losses occur through pair creation. In the 
third mechanism, plasma instability can arise from counter- 
streaming of electrons and positrons in oscillations. In the 
forth mechanism, particles accelerated in LAEW directly 
emit electromagnetic radiation, which has close analogy to 
synchrotron radiation or inverse Compton scattering. It can 
be shown that damping due to these radiation processes 
is indeed weak, with the damping time being much longer 
than the wave period and generally longer than the light- 
crossing time over the gap. Although in principle one may 
determi ne the damping from Q-dSI, averaged over a wave 
period l|Asseo. Kennel fc Pellal 19771 ). here in estimating the 
effect of wave damping, we adopt a different approach in 
which the wave damping is related to energy losses by a 
single particle. 

To estimate the damping time we consider how the to- 
tal wave energy density evolves. The total wave energy den- 
sity can be written as a sum of the electric energy den- 
sity, Ue = £qE 2 /2, and the energy density associated with 



particle oscillations, U p = m e c 2 'y m ax{n+ + ti-)Ngj. Using 
Jmax « u m ax « e 2 E 2 / (m 2 c 2 ujGj) (cf. Sec 3.2), one finds 
Up fa Ue - Let the average energy loss for a single particle be 
(7), where the average is taken over one wave period. The 
typical damping time scale can be estimated from 



Ue + U p 



2 7 „ 



(n+ + n_)A r Gj(7)m e c 2 



<7> 



(40) 



As an example, for curvature radiation, one has 7 curv = 
— (2r e c/3i?c)7 4 , where r e ~ 2.8 x 10 _15 m is the classical 
electron radius. Using the average ( 7 4 ) « 0Aj max in the 
limit Pv 3> 1, one obtains ( 7c urv) ~ —0A(2r e c/3Rc)'ymax- 
For R c w 3 x 10 5 m and ^ ma x = 10 6 , one has r ~ 
7 .5R 2 1 (r e cy max ) w 0.8 s. Thus, the damping time is much 
longer than the wave period (2-k/cj) and considerably longer 
than the light-crossing time over the gap, denoted by r g . 
In general, one has r g ^ R/c w 3 x 10 -4 s. A similar esti- 
mate can be obtained for other three processes and it can be 
shown that our assumption r ^> 2-n/u) is valid and that in 
general r > r g . For fast pulsars, curvature radiation may be- 
come efficient, as a result of a smaller curvature radius, and 
the damping time may become comparable or even shorter 
than T g , but it is still much longer than the wave period. 
Similarly, for high magnetic field pulsars (B > B c ) with hot 
polar caps, resonant inverse Compton scattering l|Sturner| 
Il995t) (cf. Sec 4.2) can be efficient and can also lead to 

T <T a . 



4 PAIR CREATION 

In this section we discuss pair creation occurring in a LAEW. 
The dominant pair production process is single photon de- 
cay in superstrong magnetic fields. We consider two main 
emission processes that produce pair-producing photons: 
curvature ra diation and resonant inv erse Compton scatter- 
ing (RICS) i|Sturnerlll995l ; lLudll996l ). The latter is inverse 
Compton scattering in cyclotron resonance, involving scat- 
tering of thermal photons from the star's surface by rel- 
ativistic electrons or positrons. Although other processess 
such as nonresonant inverse Compton scattering may also 
contribute to pair production, we only focus these processes. 
How pair creation affects the wave depends on the ratio of 
the pair production free path, \ p , and the wavelength, A. 
For Xp/X 3> 1, a photon travels many wavelengths before it 
decays into a pair. Pair injection can be treated as nonoscil- 
latory, with the injection rate derived from a sum over the 
pair production during many oscillations. In the opposite 
limit Xp/X 1, pair creation is locked in oscillations: pairs 
are injected at particular phase during each oscillation. Here 
we only discuss the first limit as it is applicable for typical 
pulsars. 



4.1 Curvature radiation 

The free path for pair production can be written as X p — 
Asi + Asp, where As; is the characteristic length for emis- 
sion of a photon at energy > 2m e c 2 and As p is the path 
length that the photon needs to travel before decay into 
a pair. For curvature radiation one can show that the for- 
mer is much shorter than the latter. To estimate Asi one 
writes the production rate, dn^ /dt, of forward (+) and 
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backward (— ) propagating curvature photons at the energy 
£ c as a ratio of the radiation power, P cur v ~ — jcurv, to the 
characteristic energy, e c , of curvature photons emitted by 
a relativistic electron moving along a curved field line with 
a curvature radius, R c . Using drv^ /dt w —jcmv/sc, where 



(3A C /2P C ) 7 3 



2Pr 



-1/2 



(7/10 6 



and A c 



h/m e c 



3.86 x 10 m is the Compton wavelength, one obtains 



dt 



R c 



(41) 



where aj w 1/137 is the fine constant. Since dn^ /dt ~ 

dn^ /dt for z R, the forward and backward components 
are approximately symmetric. For a dipole magnetic field, 
the curvature radius is given by R c = (4/3)(crP/27r) 1//2 « 

2.9 x 10 5 (r/i?o) 1/2 Po.i 2m - We have 



As, 



P c 



10 6 



(42) 



A photon with energy e 7 ~ e c > 2 needs to travel a further 
distance before being converted to a e pair. This distance 
can be estimated as follows. The opacity of a photon in a 
strong magnetic field is a function of ip = 0.5eb£ 7 sin # 7 b, 
where 9 1 b i s the propa gation angle of the photon and 
£b = B/B c l|Erber| [T966). Generally, pair creation requires 
if) ~ 1/15. A pair is produced when the opacity reaches 
unity. Using 9~,b ~ As p /R c , one has As p = 2i{)R c /eb£c 
to produce one pair. Since the maximum 7 is limited by 
radiation-reaction, denoted by jr, one may obtain a lower 
limit to the free path, 



A» 



Rc_ 

£c,R 



£c,R { 3A C 

~ \2R~ C 



1/3 



2^ 

SB 



(43) 



where £ Ci _r = (3A c /2P c )7fj. The right-hand side is ~ 500 m 
for R c = 3 x 10 5 m, ijj = 0.01, e B = 0.1, e 7 = e c , R = 10 2 . 
Using the parameters in (TjTJ one has A < A p for f3 v < 
5 x 10 2 . 



4.2 Resonant inverse Compton scattering 

Similarly, one may estimate A p for RICS. For thermal pho- 
tons at energy < I/7, one can ignore the Klein-Nishina 
effect; the pr oduction rate of the scat tered photon at energy 
e s ~ es7 is l|Sturnerlll99"5l : lLudll99d ) 



dt 



Qx^ ^Qsb ccr e ff 



87T 2 7 2 



A^ 



(44) 



where = 1.7 x 10" 4 (T S /10 6 K) is the normalized temper- 
ature T s of the polar cap, cr e g « 37T(T:r/4a/ is the effective 
cross section of RICS, and ar ' = — ln[l — exp(— £b/©7(1=f 
/3cos# m ))] with 9 m t he maximum propagation angle of 
the incoming photon (|Dermer|[l99fjh . Since ' > x^ + \ 
pair production by particles moving toward the star is 
more efficient than particle s moving away from the star 
<|Harding fc Muslimovl [i"99Sh . For both cases, the rate in- 
creases with decreasing 7. However, the process becomes 
less efficient at low energy as the cyclotron resonance con- 
dition becomes difficult to satisfy (x^ decreases exponen- 
tially when 7 is too low for the condition to be satisfied). 
This leads to an estimate of As^: 



As, w 10 e 



1/3/ T s \- 2 / 3 /z±\ 2 / 3 4/api/, 



(45) 



We write 7 as a fraction e < 1 of the maximum potential 
drop across the polar cap. Since < I/7, one has 



As P = > w ® . 

7£s es 



(46) 



For = 1.7 x 10 -4 and e B = 0.1, one has As p w 23 m. One 
concludes that for a moderate /3v > 1, A p ^> A applies to 
RICS. For 7 > 1/0, the scattering is in the Klein-Nishina 
regime, which is not discussed here. 



4.3 Pair injection 

Since for X p 2> A, one may regard dn^ /dt as a constant, 
where 7 is replaced by its average (over the period). One 
may write the source term in ([TJ as 



Q= I(a+iV + +a_iV_), 



(47) 



where a± = {dn^ /dt) + (dn^ /dt) is a constant. It should 
be noted that inclusion of higher generations of pairs is nec- 
essary in conventional treatments of pair creation in a pair 
formation front, especially for curvature radiation. Inclusion 
of higher generations would effectively increase a± in (45). 
However, the assumption that the free path for pair produc- 
tion is much longer than the wavelength remains valid. 

Assuming /3v 3> 1, integration of (|7J results in an ex- 
ponential growth 



N± ~ n± exp 



a+ + a_ 

2/V 



(48) 



where a± = a±/u>GJ- One can show that for \ ~ Pv^G-jt 
as 0v » 1, Eq (|48)l reduces to AT± oc exp[— (a+ + Q-)t/2] 
which reproduces the numerical result in iLevinson et al.l 
l|2005l ). 

We envisage pair creation as being only a minor pertur- 
bation except when the large-amplitude oscillation is being 
set up, as described in the purely temporal case by Levinson 
et al. (2005). Here one may derive the condition under which 
the effect of the source function in the equation of motion 
((8| can be ignored. This effect can be characterized by a 
parameter S± = Q/2/3vN±, where we assume (3v 3> 1. The 
pair production has only a minor effect if 8± <C 1. Using (|47[1 
and (|48[) . this condition can be writtten in the form 5± — 
(a+n+ + ct-n-) / (4f3vUGjn±) < 1. As an example, for cur- 
vature radiation one has S± ~ 2 ma ~ 1 afcj/2R c f3v^GJ <C 1, 
where one assumes there are m g generations of pairs. 



5 LOW-DENSITY LIMIT 

The low-density regime is applicable if there are initially 
insufficient charges to provide the GJ charge density. One 
can show that in this limit both monotonic and oscillatory 
solutions exist, with the former corresponding to rapid ac- 
celeration of charged particles from the surface. 



5.1 Monotonic acceleration 

A relevant example of a low density corresponds to a space- 
charge-limited flow (SCLF) from the surface, which has 
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so far been discussed only in the context of the steady 
state jArons fc Scharlemannlll979l ; iMestel fc Shibatalll994l ; 
lHarding fe Muslimovlll998l ). The case of a vacuum-like field 
is applicable for puls ars with Q, ■ B < w hen ions are tightly 
bound to t he surface jMedin fc Laill2007l ). This case was con- 
sidered in ILevinson et alJ <|2005h ~ and is not discussed here. 
We consider an outflow of electrons from the polar cap with 
fl ■ B > 0. The initial electric field at the surface is assumed 
to be small. Assuming n+ = 0, one obtains 



$ = El + 2 (7 - j )t]gj - (u - «o)io|| 



(49) 



where (|14|) is used. A solution for $ > is j i< < rjcj for 
u > uo ^ 0. As shown in figure[8] u only has a lower bound, 
corresponding to the positive part of the solid and dashed 
lines, and in these cases electrons can be accelerated mono- 
tonically outward. When j u > -qcj (dotted line), u has both 
upper and lower bounds and the solution must be oscillatory 
with the upper bound being well below the pair production 
threshold. Such lo w amplitude oscilla tory solution was also 
discussed recently l|Beloborodovll2007h . It should be pointed 
out here that the effect of the conducting wall is not con- 
sidered but can become dominant in determining the accel- 
eration (see further comments in Sec 5.2). For (3v ^> 1 one 
has d£ ~ ~Pvdu in (|22[) . which gives rise to a monotonic 
solution for the electric field: 



£11 



&o - -3-X, 
Pv 



(50) 



with Sr/ — t\gj — Jo|| > and Eo < 0. The electron's mo- 
mentum increases according to 



, \ E "\ , *V 2 
Pv lp v 



(51) 



A special case is an initial phase \ = chosen corresponding 
to t — and z = 0, i.e., at the surface where one has Eo « 0. 
We envisage that E\\ grows to reach the value at which pair 
production occurs and the oscillatory phase takes over with 
Eo in Sec 3 replaced by the threshold Eu . The acceleration 
discussed here is different from the usual SCLF models; the 
electric field is predominantly inductive. 

In the temporal case, an electric field increases linearly 
with time and thus, the 4- velocity of the particles varies 
quadratically with time: tt_ oc SrjP /2 with i = LUGjt. The 
time for an electron to be accelerated to the pair production 
threshold ^f t h is 



t = 



7th 
|J77| 



1/2 



(52) 



The time (|52|l should be limited by the travel time of a 
particle across the acceleration region. Assuming the region 
to have a size of AL, the limit gives rise to a condition 
t < AL/c for a pair cascade to occur. 

In general, the relative importance of acceleration due 
to inductive and non-inductive (potential) electric fields can 
be determined from (u_)i a ductive/(w-)static ~ Pv- F° r Pv ^ 
1, inductive acceleration is dominant. The acceleration can 
be regarded as purely electrostatic only in the special case 
Pv — > 0, which is discussed in 5.2. 




Figure 8. Plot of $ as a function of u in the low-density limit 
for (VGJ,jo\\) = (-1>-1-1) with «-0 = 2 (solid), (r]GJ,3o\\) = 
(-1,-1.1) with m_ = 0.1 (dashed) and (r)GJ,jo\\) = (-1,-0.9) 
with u—o = 0.1 (dot). The third case only permits an oscillatory 
solution with u being limited by u max ~ 8.5. The initial electric 
field is assumed to be Eq = 0. 



5.2 Comparison with the steady-state models 

The steady state can be regarded as a limit Pv — » in 
which the parallel electric field is purely static. For Pv — > 
the phase depends on spatial coordinates only and can be 
written as \ = —lugjz/c = —z. Substituting (149 \ into (|22l) 
and using d£ m dr/, one obtains 



(53) 



where Eo < 0. The monotonic solution (|53[1 can be regarded 
as an oscillatory solution in the long period limit T — > 00. 
From the current-charge invariant (|13[) one has t](0) — t)gj = 
— (P-on- + jo\\)/Pv- Since r;(0) — t\gj must remain finite, 
the limit Pv — > implies that the initial current is j\\ (0) ~ 
—P-on,-, which matches the constant current j||(0) = joy. 
The electron's momentum is derived as 

u- = u-o + \Eo\z + \Srjz 2 . (54) 

As for (I50p . for electrons to be accelerated outward one must 
have Sr/ > 0. That the same condition (Sr/ > 0) is required 
in both cases is hardly surprising. In (|50|l . one has d£ ~ 
—Pvdu < and \ > 0, while in (|50[) one has d£ ~ du > 
but x < 0. If Eo ~ 0, the electron's momentum increases 
with z quadratically, U- ~ 8r/z 2 /2. 

In the c onvent i onal SCLF model s 

l|Arons fc Scharlemannl 1 19791 : lHarding fc Muslimovl Il998l ). 
En < is obtained with 5rj < 0, by imposing an upper 
boundary, usually located at the PFF, where En = 0, and 
a conducting surface of the side wall of the open field line 
region; in these models jo\\ is then determined locally by 
these boundary conditions. The basic assumption in the 
SCLF models is the nonconstancy of Sr/ along flow. This 
means that if one sets 8rj = initially, a nonzero Sr/ 7^ 
develops along the flow inducing a parallel electric field. 
Two effects that lead to Sr/ < have been considered in 
the literature, including field line curvature, corresponding 
to the field lines curving toward the rotat ion axis, and 
frame dragging (|Muslimov fc Tsyganl Il992l ). The latter 
dominates near the star; the effective angular velocity, so is 
the GJ density, is reduced by a factor (1 — k g (R/r) 3 ) < 1 
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as compared to that observed in a flat space at infinity, 



where k g = 2GI/{c I R i 



0.15/38, z < R 



It) 4 



and 



7 3 s = 7/(10 JM kg m ) is th e moment of inertia of the star 
|Muslimov fc Tsygarj[l992T ). If one assumes Sr] = initially 
at the surface, one has Srj = 3k g r]Gjz/R < for r/cj < 0. 

Eq l|53|l a nd (|54[) are similar to the result derived by 
IShibatal |l997h based on a generic SCLF model in which 
no specific local boundary condition is imposed. When j u 
is treated as a free parameter, for initially 8r/ — 0, 8-q > 
is required to produce E\\ < 0. This can occur only on the 
curving-away (fr om the rotation axis) field l ines along which 
decreases (|Shibatall 19971 ; lMestel|[l999l ) . A major prob- 
lem with this scenario in the context of the stea dy-state 
limit is that the g rowth in \E\\\ is unstoppable (jShibatal 
1 19971 ; iMest cl 1999). However, such run-away growth does 
not occur in our oscillatory model because pair creation ul- 
timately leads the system to switch to an oscillatory phase, 
as discussed in Sec. 3. For Sr] < 0, Eq (|54p implies an os- 
cillatory (in space) solut i on similar to t hat found previously 
jMestel fc Shibatalll994l ; IShibata|[l997l ). When the acceler- 
ation region extends to > R(R/ Rlc) 1 ^ 2 , the effect of the 
conducting side wall, at which En = 0, becomes important. 
When such effect is included acceleration of outflowing elec- 
trons is possible at > R(R/ Rlc) 1 ^ 2 even when Sr/ < (pro- 
vided that an electric fie ld arising from such effect dominates 
over that from Sr] < 0) (| Shibat al 1 1 997l ^ . 



6 CONCLUSIONS AND DISCUSSION 

We present an oscillatory polar gap model, in which the sys- 
tem initially undergoes a low-density phase, involving rapid 
acceleration of particles to ultra high energy, initiating a pair 
cascade. The system evolves to an oscillatory phase. The 
oscillations are treated as a superluminal, large amplitude 
electrostatic wave that propagates along the magnetic field. 
The charge continuity equation implies a current-charge in- 
variant (j|| — f3vrj — const) that is independent of pair cre- 
ation. As a result, the phase velocity f3y is no longer a free 
parameter and can be written in terms of the initial velocity 
and density of the plasma. It is shown that only the superlu- 
minal case (5v > 1 is relevant here. An analytical formalism 
for LEAWs is derived in the high-density regime in which 
the pair density is higher than the GJ density. We ignore 
wave damping in our analytical solution. Neglecting damp- 
ing is justified as the typical damping time due to energy 
losses through radiation is much longer than the wave pe- 
riod. In most cases, the damping time is also longer than 
the light- crossing time over the gap. 

The model predicts an outflow of relativistic pairs due 
to particles being dragged along in LAEW. Such feature is 
needed to avoid overheating of the polar cap. Outflowing 
pairs would contribute to the pulsar wind. Pairs oscillate 
with a net drift velocity directed along the magnetic field, 
producing a current that oscillates about the global constant 
current jo . The amplitude of the oscillating current is larger 
than the global current by a large factor that is of order of 
magnitude the ratio of the pair density to the GJ density. 
The wave form of an inductive electric field is characterized 
by a triangular shape, which can be understood as the cur- 
rent being nearly constant except for a brief period during 
which it switches sign. The basic features of the oscillations 



are not sensitive to the initial conditions including the elec- 
tron's or positron's initial velocity. 

There are two possiblities for particle acceleration in 
the initial phase that leads to oscillations: (1) a vacuum-like 
initial electric field, which may be applicable for the polar 
cap where charges are tightly bound to the surface, and (2) 
SCLF, in which there is a n ample supply of ch arges. The 
first case was discussed in iLevinson et all l|2005h . Here we 
consider specifically the SCLF case where an initial elec- 
tric field appears as a result of an imbalance between the 
charge density and the GJ density with the latter mimick- 
ing the positive background charges. Electrons are acceler- 
ated monotonically in the electric field that increases lin- 
early with the phase x- Since \ comprises both temporal 
and spatial variables, such particle acceleration arises from 
a mixture of inductive and non-inductive effects. An inter- 
esting limit is V — > oo, in which the electric field becomes 
purely inductive. Qualitatively, the usual steady-state the- 
ory can be reproduced in the limit of a zero phase speed. 
In this limit, the system is time independent and the ac- 
celeration occurs at a specific spatial location. By contrast, 
acceleration due to an inductive field can occur everywhere 
in the region concerned. 

An implication of the oscillatory model is the predic- 
tion of plasma instability arising from counterstreaming of 
electrons and positrons; in each oscillation electrons and 
positrons are accelerated in opposite direction and such 
counterstreaming provides an ideal condition for two-stream 
instabilit y which may be directly relevant for pulsar radio 
emission l|Verdon fc Melrose 2007). Although various forms 
of streaming instability have been discussed in connection 
with the radio emission in conventional models, the growth 
rate is generally too low to be effective, requiring some sep- 
arate assumption to enhance it. In the oscillatory model, 
the relative streaming of electrons and positrons allows the 
maximum possible growth rate for the two-stream instabil- 
ity, at phases where the counterstreaming is nonrelativistic 
or mildly relativistic. Apart from the two-stream instability, 
LAE may also operate in conversion of LAEW to electro- 
magnetic radiation. 

In an oscillatory pulsar magnetosphere, cyclotron res- 
onance can have a significant effect on the propagation of 
the coherent radio emission. In the conventional polar cap 
models, for a wave propagating at an angle 9 to the mag- 
netic field, the cyclotron resonance occurs preferentially in 
the large-angle regime 6 ^> 1/7 at a frequency oj — Qe/jO 2 , 
located at a radius, which is generally in the outer magne- 
tosphere, r c ~ (Q e o/ , yQ 2 u j) 1 / 3 , where O e o is th e cyclotron 
frequency at the surface (|Luo fc M elrose 2001). In the os- 
cillatory pulsar magnetosphere, the cyclotron resonance can 
occur at uo = £l e /2'y for particles moving toward the star. 
The cyclotron radius r c varies with oscillating 7, with the 
smallest radius being (r c ) m j n w (f2 e o/27max^) 1 ' /3 ; This ra- 
dius is smaller than in the usual polar cap models by a factor 
(6» 2 /2) 1/3 « 0.17 for 9 = 0.1. 

There are some limitations of our model, notably the 
one-dimensional assumption that may not be realistic for 
an acceleration region extended beyond > R(R/ Rlc) 1 ^ 2 ■ 
The effect of the side wall of the open field line region needs 
to be included in the calculation. Such region can be mod- 
eled as a wave guide and propagation of LAEWs in such 
wave guide will be discussed elsewhere. Nonetheless, from 
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this one-dimensional, analytical model we are able to derive 
some fundamental features of LAEWs that should remain 
valid qualitatively for a more general, three-dimensional case 
as well. The fluid treatment adopted here may not be accu- 
rate for pulsar plasma as numerical simulations showed that 
pairs from a cascade generally have a broad distribution 
llArendt fc Eilekl 12002) and inclusion of a particle distrbu- 
tion requires a kinetic formalism which is beyond the scope 
of this paper. 
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APPENDIX A: DERIVATION IN THE 
TEMPORAL GAUGE 

In this appendix, we outline an alternative derivation of the 
wave equation in the temporal gauge. Electric and magnetic 
fields can be expressed in terms of a vector potential 



E 



M B 

dt 

From ^ we have 



V x A. 



V(V.A)-V 2 A + 1^ 



Mo( J - 



(AI) 



(A2) 



Assuming that all the relevant quantities are functions of X: 
the parallel (to k) component of (|A2|) takes the form 

(A3) 



with An = K ■ A and 

Joy = V||(V • A) — V 2 A|| 

= Vi|(Vx • A x )- VIAn 



(A4) 



The equation of motion ((2} can be written into a similar 
form to (fT7)> : 



spv 



+ 



11 dx m e c dx 

which can be integrated to yield 

l-pvps , _ se/3 v 



C UeC 2 2iV s 



,(A5) 



(1 - A 2 ) 1 / 2 



(A|-^o||) 



-- / X f^2 ~ Sr u °) d x, (A6) 

c J Q \m e c 2 2N 3 J 



where A ^ = A^(x = 0) and s — ± corresponds to electrons 
(+) and positrons (— ). Assuming A — (e/m e c)(A\\ — Aoy), 
for the electron component /3 = j3~, one obtains 



dA 



(1-/32)1/2 



n_+(??? + -(l-/3v/3) 



Joy 



(A7) 



where / is an integration constant. Clearly, a physical solu- 
tion requires the RHS to be non-negative; this condition can 
be satisfied only if /3 min < f3 < /3 ma2 , , where f3 min and /3 mol( , 
are the minimum and the maximum velocities at which the 
RHS is zero. The electric field can be found from 



-Pv 



dA 



(A8) 



Eq (1 2 2 tl and (1231) r eproduce an an alytical form similar to 
that given bv lAkhiezer et all (|l975h : 



(1-/3 2 .) 1 / 2 (l-/? 2 ) 1 ^ 



-1/2 



de = 2 i/^ X; (B2 ) 



E» = ±2 



1/2 met 



(1-/3 2 ,) 1 / 2 (1-/3 2 ) 1 / 2 



1/2 



(B3) 



Consider the relativistic limit 7 m = 1/(1 — Pm) 1 ^ 2 3> 1; 
Using /3 « ±(1 - 7~ 2 /2), one finds 



1 ± /3v 1 + 



7m 



+ 



2 T Pv 1 + 



(7m-l) 1/2 

(7- m -7) 1/2 



Pv 



arctanh I 1 — ] _ arctanh f 1 



7m 



7m 



X- 



(B4) 



Since /3 is a periodic function of x, we can define a period T 
by 



0v -P 



1 



-1/2 



Jl-Z^m) 1 / 2 (l-/3 2 )l/2_ 

= V2u) p p v T. (B5) 
The RHS can be written into the form in the 7 m 3> 1 limit: 



Pdj 



87 



1/2 



/x (7m -7) 1/2 

If we define a frequency u = 2n/T, (|B5|) leads to 
\/27rWp 



1/2 



47 

which is similar to J33 



(B6) 



(B7) 



APPENDIX B: ELECTRON GAS 

Here we reproduce the known result for a LAEW in an elec- 
tron gas, by retaining the electron component only. Specif- 
ically, we set Po = P~o = and 77+ = 770 = r\ G] = 0. The 
maximum velocity can be expressed in terms of Eq : 



